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Abstract 

We reconsider the tree level color-singlet contribution for the inclusive J/ ip production in T decay 
with the ttj order QCD process T — > J/t/j+cc+g and o?o? s order QED processes T — > 7* — > J/tp+cc 
and T — > J/tp + gg- It is found that the contribution of the QED process is compatible with 
that of the QCD process, and the numerical results for the QCD process alone is an order of 
magnitude smaller than the previous theoretical predictions, and our theoretical prediction in total 
is about an order of magnitude smaller than the recent CLEO measurement on the branching 
fraction B(T — ► J/tp + X). It indicates that the J/tp production mechanism in T decay is not well 
understood, and further theoretical work and experimental analysis are still necessary. 

PACS numbers: 12.38.Bx, 12.39. Jh, 13.20.Gd 
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I. INTRODUCTION 



Since the discovery of cc state J/ip and bb state T more than three decades ago, heavy- 
quarkonium system has served as a good laboratory for testing QCD from both perturbative 
and non-perturbative aspects. With the accumulation of new experimental data and the 
development of interesting theory, considerable attention has been attracted to study heavy- 
quarkonium spectrum, decay and production (for a review see 1|). 

On the theoretical side, the non-relativistic QCD(NRQCD) 2] effective field theory was 
introduced, based on which the production and decay of heavy quarkonium can be calculated 
with a rigorous factorization formalism. This formalism separates the physics on the energy 
scale larger than the quark mass mg, related to the annihilation or production of QQ pair, 
from the physics on the scale of rriQV 2 order, relevant to the formation of the bound state. 
Consequently, the inclusive production and decay rates of heavy quarkonium are factorized 
into the product of short- distance coefficients, which could be calculated perturbatively 
as the expansion of a s , and the corresponding long-distance matrix elements, which are 
determined by some non-perturbative methods. The long-distance matrix elements are 
weighted by the powers of v, the velocity of heavy quark in the rest frame of the bound 
state. One important feature of NRQCD is that it allows the contribution of QQ pair 
in color-octet configuration in short distance, and the color-octet state will subsequently 
evaluate into physics state through the emission of soft gluons. 

The introduction of NRQCD has greatly improved our understanding of the production 
mechanism of heavy quarkonium. One remarkable success of NRQCD is that the transverse 
momentum (p t ) distributions of J/tp and ip' production at Fermilab Tevatron[3] could be 
well described by the color-octet mechanism|4]. However, this mechanism could not correctly 
explain the CDF measurements of J/ip polarization 5|. Just about one or two years ago, 
the next-to-leading order (NLO) QCD corrections to both the color-singlet and color-octet 
processes have been obtained. For the color-octet process^, it is found that the leading 
order (LO) results are little changed when the NLO QCD corrections are taken into account. 

In the color-singlet case, the theoretical predictions at QCD NLO are significantly changed 

n 

from the LO results on the pt distribution and polarization of J/^[7|. Although this still 
could not resolve the puzzle. The large impact of the color-singlet NLO QCD corrections 
on the LO results indicates that the contribution of the color-octet mechanism may not 
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as important as we expected before. Furthermore, the theoretical predictions [8(1 for the 
Pt distribution of T can compatible with the data of T production at Tevatron|9j within 
the theoretical uncertainty when considering some of the important next-to-next-to-leading- 
order (NNLO) aP s contribution. However, it still cannot explain the recent polarization 
measurement by DO Collaboration [lOj 

In the case of J /if) production in e + e~ annihilation, the existence of color-octet mechanism 
also faces to a challenge. The NRQCD approach predicts that the J /if) production in e + e~ 
annihilation at LO in a s is dominated by e + e~ — > J /if? + gg, and e + e _ — > J /if) + cc and 
e + e _ — > J /if) + g, in which the first two are color-singlet subprocesses and the last one 
is color-octet subprocess. The color-octet contribution ll] predicts there is a peak in J / if) 
momentum spectrum near the kinematic end point. Unfortunately, The peak was not found 
in the experimental observation of BABAR[l2| and Belle 13|. By using the soft-collinear 



NLO QCD calculations shown that both cr[e + e 

j/^+x non _ cc -]Q,[i? 



effective theory (SCET), the color-octet predictions [14| could be softened, but it depends on 
a unknown non-perturbative shape function. Belle also extended their analysis by deriving 
associated J /if) production with cc pair from inclusive J /if) production production 151 ] . The 

> J /if) + cc + X]\m, |l3] and a[e+e~ -> 
, ll9j may be explained by considering only the contribution of color-singlet 
process. However, it point out in Ref. [17( that the color-octet contribution is still not yet 
completely ruled out due to the incomplete measurement in the experimental analysis. 

To improve our understanding of J /if) production mechanism, it was proposed 2(], 21] that 
the T decay may provide an alternate probe of J /if) production in rich gluon environment. 
Experimentally, the branching ratio of T — > J/if> + X has already been reported to be 
(1.1 ± 0.4 ± 0.2) x 10~ 3 by CLEO based on about 20 events in Ref.[22|. The ARGUS 



Collaboration obtained an upper limit of 0.68 x 10 -3 23j at 90% confidence level. With about 
35 times larger data sample than previous work, an improved measurement of J /if) branching 
ratio and momentum spectrum have been obtained recently by CLEO Collaboration with 
B(T -> J / if) + X) = (6.4 ± 0.4 ± 0.6) x 10- 4 [24]. Theoretically, the color-octet prediction 
is B(T -> J / if) + X) ~ 6.2 x 10- 4 [21] with 10% contribution from ^(25) feed-down and 
another 10% from Xc R However, it ™i found that the branching ratio of color-singlet 
process T — > J/ip + ccg is about 5.9 x 10 4 ]26], which is also in agreement with experimental 
measurement. Although both the color-singlet and color-octet decay modes may explain 
the total decay rate independently, their predictions on the J /if) momentum spectrum are 
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significantly different. The maximum value of J/ip momentum in the color-singlet and 
color-octet process are 3.7 GeV and 4.5 GeV respectively. The CLEO collaboration found 
that the experimental result of J/ip momentum spectrum is much softer than color-octet 
predictions and somewhat softer than color-singlet predictions. The process T — > J/ip + X 
also was studied in color evaporation model 27| more than thirty years ago, but this model 
can not give systematic predictions of J/ip production. Another early theoretical work on 
the process T — > J/ip + X could be found in Ref . [3] . 



There is a very we 
mental measurements 



1 agreement between the LO color-singlet predictions 26] and experi- 
24}]. But it seems difficult to understand the situation in comparison 
with the case of the J/ip production at B factories, where there are huge discrepancies be- 
tween the LO theoretical predictions and the experimental measurements. Therefore, we 
re-calculate the branching ratio of color-singlet process T — ► J/ip + cc + g in this paper. 
And the results show that it is an order of magnitude smaller than the previous theoretical 
prediction 



26]. Therefore, there is an order of magnitude discrepancy between the LO 
theoretical prediction and experimental measurement for T — > J/ip + X now. To further 
clarify the situation, we also estimate the leading-order contribution of the QED processes 
T — > 7* — ► J/ip + cc and T — > J/ip + gg at a 2 a 2 s order, in which the process T — > J/ip + gg 
includes two gauge invariant subsets, T ^7* — > J/ip + gg and T — > ^*gg followed by 
7* — > J/ip. The final results show that the contribution from the QED processes are com- 
patible with that from the QCD process. 

The rest of paper is organized as follows: In section II, the basic formula and method used 
in the calculation are presented. In section III, we describe the calculation on the branching 
ratio of the QCD process T—>J/ip + cc + g and J/ip momentum spectrum. In section IV, 
we estimate the contribution of the two QED processes T — ► J/ip + cc and T — > J/ip + gg. 
The final results and summary are given in the last section. 



II. DESCRIPTION OF OUR BASIC CALCULATION FORMULA 

At leading order in vq, for S-wave heavy- quarkonium production and decay, the color- 
singlet model predictions are equal to that based on NRQCD effective theory. Then we 
express dT(T — ► J/ip + X) as: 

dr(T ^ J/iP + X) =dt(bb[ 3 S 1 ,l] cc[ 3 S 1 ,l}+ X)(T\0 1 ( 3 S 1 )\T)(Ot{ 3 S 1 )}, (1) 
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where dT(bb[ 3 Si, 1] — > cc[ 3 Si,T\ + X) represents color-singlet bb pair in spin-triplet state 
decay into color-singlet cc pair in spin-triplet state with anything, which is calculated per- 
turbatively, and (T|0 1 ( 3 S 1 )|T) and {Of ( 3 ^)) are the long-distance matrix elements, which 
can be related to the nonrelativistic wave functions as: 

(T\0 1 ( 3 S 1 )\T) ~ ^\Rr(0)\ 2 , (Ot( 3 S l )} = |^(0)| 2 . (2) 



We employ spinor projection method 29j] to calculate the short-distance part dT. In the 



nonrelativistic limit, the amplitude of bb[ 3 Si,l\ — > cc[ 3 Si, 1] + X could be written as[30|: 



MibbfS!, l](po) - cc( 3 S, 1)^) + X) = E E E 
x ( Sl ;s 2 | lS z >(3i;3j | 1) x (s 3 ;s A \ lS z )(3k;3l | 1) 

x M(h(^, s 1 )b j (^, s 2 ) -> c fc (|, S3 )Q(y , + X) (3) 

where (3z; 3j | 1) = 5 i:j /y/W c , (3k; 31 | 1) = 8 kl /^/N~ c , (s 1 ;s 2 \ IS Z ), and (s 3 ; s 4 | 1S" Z ) are 
the SU(3)-color, SU(2)-spin and angular momentum Clebsch-Gordan (C-G) coefficients for 
<5Q projecting on certain appropriate configurations at short distance. At leading order 
in vq(Q = b, c), the projection of spinors u(^-, Si)v(^-, s 2 ) and v(^-, s 3 )u(^, S4) could be 
expressed as: 

S1,S2 

n c = E( s i; s * I i4Mj, ^(|, , 4 ) = ^fiftsjfa + 2m c ), (4b) 

where e(S" 2 ) is the polarization vector of the heavy quarkonium. For a spin=l state with 
momentum p, the sum over its all possible states S z is 

^ a \o z )tp\o z ) = y-y al 3-t- 

s z 

According to the spinor projection method, the relation between dT and \Ai\ 2 for the 
bb[ 3 S u l] -> cc[ 3 S u l}+X is 

1 1 V I A// 1 2 

rfrws,, i] - ccf 5l , i] + cc + 9 ) = -— g-MJ^ rf.. ( 6 ) 

where means to sum over all possible polarization states of the particles in this process 
and $ n is the n-body phase space. The factor (l/2iV c ) 2 with N c = 3 comes from the 
normalization factor of the NRQCD 4-Fermion operator. 



£^W) = (-&* + ^) ( 5 ) 



Since our calculation gives different results from the previous theoretical prediction [261 ]. 
we further checked our results by using two different way to do all the calculations. One 
is to apply the above formula to write a piece of program to do the calculations for each 
process described in the following two section. Another is just using the Feynman Diagram 



Calculation (FDC) Package [3l|] to generate all the needed Fortran source and then do the 
numerical calculation. We obtained exactly the same results by using these two methods. 
Moreover, to check gauge invariance, in the expression of FDC version source, the gluon 
polarization vector is explicit kept and then is replaced by its 4-momentum in the final 
numerical calculation. Definitely the result must be zero and our results confirm it. 

III. THE QCD PROCESS T — > J/ip + cc + g 

Now we proceed to calculate the total decay rate of T — > J/ip + cc + g and its contribution 
to the J/ip momentum spectrum. At leading order in a s , there are six Feynman diagrams 
which are shown in Fig. 1. The amplitude M. could be factorized as: 

M(bb[ 3 S u l](p ) -> cc[ 3 S h l](pj + c(p 2 )c(p 3 ) +g(p A )) = 

MbibbfSx, 1] - g*g*g) x Me(ffV - cc[ 3 S 1 , 1] + cc), (7) 

in which the later one is universal for all the six diagrams and it is 

M ° = ( M /9W M /,^ W»(ft). (8) 

(pa +Pi/2y{p 3 +pi/2) 2 
The amplitude of Ji4b(bb[ 3 Si, 1] — > g*g*g), for example for the first diagram, is 



' * ' 1 ' 7 'W2 + p,/2+rf-^ 7 W2-P.) ! -^' ,) 



where C\ is the corresponding color coefficient and ^ 3 is the polarization vector of the real 
gluon. The amplitude hA\ for the other five diagrams could be obtained in a similar way. 
An analytical expression of ^2 \M\ 2 is obtained in the calculation, but is too lengthy to be 
presented here. 

The four-body phase space $4 for bb[ 3 Si, 1] — > cc[ 3 Si, 1] + cc + g is defined as 

rf$ 4 (p - Pl +p 2 +p 3 +pj = n 7^i^( 27r ) 4(54 ^ - Eft) ( 10 ) 

fc=l ^ ' fe fc=l 
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FIG. 1: The six Feynman diagrams for the short-distance process: bb[ 3 Si, 1] — ► cc[ 3 5i, 1] + cc + 5. 



There are many ways to perform the four-body phase-space integration. Here we briefly 
introduce our methods. Using the two following identical equation 



d 4 P^ 
(2tt)< 



(2tt) 5 (p 234 -p 2 -p 3 -p 4 ) = 1, 



d*P. 



14 (2*)W(p M - p 3 - p 4 ) = 1, (11) 



(27T) 



we transform the four-body space into the combination of three two-body phase spaces, 
which is given by 



d$ 4 (p ~>Pi +P 2 +P 3 +P4 



iJ 2tt 2tt 



(12) 



where s 234 = P 2 34 ,s 34 = p 2 34 . The three two-body phase spaces integration are described by 
the three-momenta p 1 ip^ : p^ f * and their solid angle element dfl , d£l* , dfl** in the rest frames 
of p , p 234 , and p 34 respectively. Then the expression of four-body phase space becomes 



2tt J 8(2n) 2 m b J 2vr ' 234 



IK 



4(27t) 2 ^7 4(2tt) 



IpTI 



: dn«. (13) 



where [pj|, |p~^| and |p^*| are given in the equations below in the rest frame of p , p 234 and 
p 34 respectively 

« /ifim^ -I- ( — 4m 2 —I— jej "l 2 — Rm 2 f 4m 2 -I- .« ^ 

(14a) 



|pi| 



yi6m| + 


{-Ami + s 234 ) 2 - 


8m2(4m2 + s 234 ) 


4m& 


Y ( S 234 ~~ ( 


m c y/^TT) 2 ) ( S 234 


- ( m c + v^T) 2 ) 










2 

1^**1 _ S 34 _ m c 

34 





(14b) 
(14c) 
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TABLE I: The values of f(r) for different r = m c /rrib 



r 


0.275 


0.296 


0.317 


0.327 


0.338 


0.361 


0.381 


f(r) 


0.904 


0.567 


0.345 


0.269 


0.202 


0.105 


0.055 



The integration ranges of s 234 and s 34 are 

Am 2 c < s 234 < ( 2m b ~ 2m c ) 2 ,m 2 c < s 34 < {^/F 2 



(15) 



For space-symmetry, dQ and d<p* 34 could be integrated out directly then \A4\ 2 only 
dependent on five variables s 234 , s 34 , 9* 34 , 9**, and 0**. To get the total decay rate, the 
non-trivial integral with these five variables is performed by three steps. First, we do the 
integration dQ** in the rest frame of p 34 , then we integrate out s 34 and 8* 3i in the rest frame 
of p 234 , the last variable s 234 is integrated out in T rest frame. Since \pi\ only depend on s 234 , 
the J/ip momentum spectrum could be easily obtained by replacing ds 234 with ^#y<i|pi|. 
The phase space integrations for the total rate and J/ ip momentum spectrum are calculated 
numerically. 

By dimension analysis, it is easy to represent the decay width and differential decay width 
of T — »■ J/ip + cc + g as 

(T\0 1 ( 3 S 1 )\T) (Oti^)) 



T(T -> J/^ + cc + g) 



a:' 



ml 



■fir) 



dY 



(T -> J/^j + cc + 



a 



\g{r, \pi\/m b 



2N C 3x2N c 
(T\0 1 ( 3 S 1 )\T) (Ot( 3 S 1 )) 



(16a) 
(16b) 



where r = m c /mt, and f(r) are dimensionless, and f(r) function is same as h(r) in Ref. 26] 
To ensure the validity of our calculations, we use two different kinds of computer codes for 
cross check and obtain exactly the same results for f(r) and g(r, |pi|/m&). When r = 0.327, 
the decay width is 

a 5 s (T|0i( s Sl)|T) (Of( 3 S 1 )) 



r(T -> j/^ + cc + 



Ml 



2A^ 



3 x 2N C 



x 0.269. 



(17) 



To compare our results with those in Ref. 26j, the numerical results of f(r) in the range 
of 0.275 < r < 0.381 are listed in Tab. [I]. It is easy to see that the results of f(r) are 

n 

about an order of magnitude smaller than that given in Ref. [261] and f{r) changes a little 
sharper than that when r goes from 0.275 to 0.381. Besides f(r), the decay width T(T — > 
J/ip + cc + g) is also dependent on the choice of the values of the two long-distance matrix 
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elements (T|d?i( 3 <Si)|T), (Of ( 3 Si)), the coupling constant a s and the mass of b-quark. To 
reduce the uncertainty of theoretical predictions, we normalize it by the decay width of 
T — > light hadron, which includes two dominate decay modes T — > ggg and T —>■ 7* — ► 
qq (q = u, d, s, c). At leading order in a s and v b , we have 

20a 3 (rr 2 -9) 



F(T - ggg) 



r(T 



2A3m 2 b 



-XTIOiC^OIT), 



2irN c e 2 a e 2 h a 2 , , ^ 

9 6 -(T|0 1 ( 3 S' 1 )|T). 



Then the normalized width r^ r is given by 



-pccg 
1 Nor 



3(2AT C )2( J|a3 (7r2 - 9) + V 2^^^^ 



r/r, 



and the branching ratio turns to be 



(18a) 
(18b) 



(19) 



B(T -> J/^ + cc+(?) 



r£ x B(T 



light hadron). 



(20) 



Since the process T — > J/ip + cc + g can be viewed as T — > gg*g* followed by 
J/ip + cc, as suggested in Ref . [20| . it is reasonable to chose a s {2m c ) = 0.259. Using e u = 

0.327, m b = 4.73GeV, \RM\ 



ed = - 



3 ,e s 



■i e 

3' Cc 



3 , e b 



3' 



128' 



1.548 
4.73 



0.81GeV 3 being calculated in potential model 32] and B(T — > light hadron) = 92% 34], we 
predict 

i3(T -> J/^ + cc + ^) = 2.12 x 10" 5 (21) 

The normalized J /if) momentum spectrum dY 'not/ 'd\pi\ is shown in Fig. 3. It is easy to see 
that the shape of the J /if) momentum spectrum is similar with that in Ref. 26|], although the 
prediction for the total decay width is an order of magnitude smaller than the experimental 
data. 



IV. THE QED PROCESS T — > J /if) + X 

There are two QED processes T — > J /if) + cc and T — > J /if) + gg at the leading order in 
a s and a. Both of them are considered in this work. We will present a few simple steps and 
analytic results for them in the following. 



9 




WWWV^^Jj 

(cl) 



ccpSi, 1] 



3 S!,1] 




(c2) 



FIG. 2: The typical Feynman diagrams for the QED processes of inclusive J /ip production: (a) 
bb[ 3 S l , 1] - 7 * - cf Si , 1] + cc, (b) bb[ 3 S l , 1] - 7 * - c[ 3 Si , 1] + 59, (c) 66[ 3 Si , 1] - ccfSx , 1] + <?<?. 



A. T -» 7* -> J/^ + cc 

At the leading order, there are four Feynman diagrams for T(po) ~~ > 7* ~~ > ■VV'CPi) + 
c(p2)c(ps), two of which are shown in Fig. 2a. The calculation procedure for this process is 
very similar to that for the production in association with cc pair in e + e~ annihilation. 
The differential decay width is given by 



dY 

d \Pi\ 



(T -> 7 * -> J/^ + cc) 



27rC A ^cgc^ 2 ^(T|0 1 ( 3 g 1 )|T)((9f( 3 g 1 ))v / ^-4^ 
9(2A c ) 2 m6 r xi 4 (k - Xi) 3 (-2 + x x ) 2 (k + Xi) 3 



(2 rexi(-2re 6 (l + 2r 2 )xi 2 + /t 4 (6r 6 (-4 + 3xi 2 ) + 2x x 2 (-4 + x x 2 (-2 + 9xi)) - 

4r 4 (16 + xi(-16 + Xi(-8 + 9xi))) + r 2 (-2 + xi)(16 + Xi(-24 + x x {-U + 39xi)))) + 

2k 2 xi 2 (8xi 2 + 7xi 4 - 18xi 5 - 4r 6 (-8 + xi(8 + x x )) + 4r 4 (20 + x^-40 + a* (13 + 4xi))) + 

r 2 (32 + xi(-96 + a* (60 + (76 - 37x x )xi)))) + xi 4 (6r 6 (4 + xi 2 ) + 2xi 2 (-4 + xi 2 (-4 + 9xi)) + 

4r 4 (8 + a* (32 + (-26 + xi)xi)) + r 2 (-32 + a* (128 + x x (-124 - 60xi + 39xi 2 ))))) + 

(k - xi) 3 (k + xi) 3 (-6r 6 (4 + xi 2 ) + 2x x 2 (4 + xi 2 (-13 + 8xi))4r 4 (-16 + x x + 



(16 + Xi(-4 + 5xi))) + r 2 (-32 + x x (64 + x x {A + (4 - 7xi)xi)))) log 



X\ + k' 



(22) 
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where Ca = 3 and Cf = § are the color factors, and there are x\ = a/Ip*l | 2 + 4m 2 /mb and 
k = y/( Xl + 2r)(x 1 - 2r)(l + r 2 - rn)(l - Xx)/(1 + r 2 - xi). 

Integrating |pj numerically and normalizing T(T — > 7* — > J/ip + cc) by T(T — > 
light hadron), we obtain 

r(T -> 7*^ J/jj + cc) 3.85a 2 ag(Q?( 3 g 1 )) 
Normal r(T ^ Hght hadron) " 6iVc ( |o_ a 3 (7r 2 - 9) + £ 9 27re 2 e 2 a 2 )m 3 , ' 1 J 

By choosing the same numerical values for r, m&, e g , a (C7 ( 3 Si)) and i3(T — > light hadron) 
as those in Sec. Ill, the numerical result is 

B(T -> 7* -> J/V> + cc) = 1.06 x 10" 6 , (24) 

and the normalized J/ip momentum spectrum is shown in Fig. 3. 

B. T^J/^ + gg 

The process T(p ) — > J/^^J + g(p 2 )g(p 3 ) includes two parts, T — > 7* — > J/ip + gg and 
T — > (?(?7* and 7* — >■ J/ - ?/'- There are six Feynman diagrams for each part at the leading 
order with the typical ones shown in Fig. 2b and 2c. To calculate the contribution of the 
two parts together, the differential decay width is represented as 

— (T - J/V> + gg) = 



dlpj 9(2Ay 2 m£r 3 X! (-1+r 2 ) (2r 2 -Xi) 3 (-2 + x 1 ) c 

((-1 + r) (1 + r) (2 r 2 - 27) (-2 + 27) ^-4r 2 + Xl 2 (8 + 8 r 8 - 4 r 6 (-4 + 3 Xi) + 

r 4 (-2 + xi) (-16 + 7m) + 27 (-12 + (7 - 2m) 27) - 2r 2 (-1 + 27) (8 + (-7 + 27) 27)) + 

2 (1 + r 2 - 27) (-((2 r 2 - 27) (8 + 2 r 8 + 27 (-12 + 5 27) + r 6 (40 + 27 (-32 + 5 27)) + 

r 4 (6 - (-2 + 27) 27 (-19 + 6 27)) + r 2 27 (-6 + 27 (13 + 2 (-5 + 27) 27))) 

log( ~^ + Xl ~ V ~^ + Xl J ) + r 2 (-2 + 07) (8 r 10 - 12 r 8 27 + Xl 2 (5 + 2 (-3 + 27) Xl ) + 
— 2 + 27 7- v — 4r 2 + 27 2 

r 6 (6 + m (-6 + 5 xi)) + r 4 (40 + 27 (-38 + 13 27)) + r 2 (2 + 27 (-32 + (31 - 10 27) 27))) 

-2r 2 + 27 + V-4r 2 + X1 2 

lQ g( 2j _ / , 2 ^ 2 )))> ( 25 ) 

—2 + Xi — v — 4 r 2 + X1 2 

Where there is Xi = a/|]7 | 2 + 4m 2 /m&. And the normalized decay width becomes 

r99 r(T->J/V> + gg) 60.8aW s (Of ( 3 ^)) 

Normal r(T ^ light hadron) GiVed^a 3 ^ _ 9) + ^ 2 7Te 2 e 2 a 2 )m 3 " 1 J 
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By using the same parameters as above. We obtain 

B(T -> J + gg) = 1.67 x 1(T 5 



(27) 



and the normalized J /if) momentum spectrum is plotted in Fig. 3. In the numerical result, 
about 85.2% contribution comes from the T — > gg^f*(J/if>) part, 18.2% from the T — > 7* — > 
J/ipgg part and —3.4% from the interference part. 

V. SUMMARY AND DISCUSSION 

To sum up all the contributions of the color-singlet QED and QCD processes considered 
above, the branching ratio of direct J/ if) production in T decay is 

£Dircct(T -> J /if) + X) = 3.9 x 10~ 5 , (28) 

and the corresponding normalized J /if) momentum distribution is given by the solid line in 
Fig. 3. It can be seen in Fig. 3 that the contribution of the QCD process is dominated in 
small p,p region, while the effect of the QED process J /if) + gg is more important in large 
p<p region. In Eq. (I25I) and the dot-dashed line in Fig. 3, the logarithmic divergence at the 
kinematic end poi nt is obvious shown for the QED process J/if> + gg. It was pointed out 
in Ref Q, 0] that both the a s and Vb expansion failed near the kinematic end point 
region in the similar processes e + e~ — > J /if) + X and T — > 7 + X because of the large 
perturbative and non-perturbative corrections, and the logarithmic divergent behavior can 
be soften by applying the resummation in the SCET. Whatever it can improve the J /if) 
momentum spectrum largely near the kinematic end point, but the corrections to the total 
decay width is small. Therefore we omit the resummation effect here. 

Our calculations show that at the leading order in a s , Vb and v c , the QCD process 
T — ► J/ip+cc+g only accounts for 54.4% of the LO theoretical prediction for total branching 
ratio, in spite of a enhancement factor a^/a 2 that is associated with the QCD and QED 
coupling constants when compared to the QED processes. The main reason lies on the 
fact that the virtuality of the two virtual gluons are both of order in the QCD process 
while the virtuality of the photon is fixed to Am 2 c in the QED processes dominated by 
T — > ggj*(J/if>), and moreover the four-body phase space of the QCD process is also less 
than the three-body one of the QED processes. 
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Total 




p(GeV) 



FIG. 3: The contributions of QCD process T — > J/ijj + cc + ^(dashed line) and QED processes 
T — > J/iJj+gg (dot-dashed line) and 5 times of T — > 7* — > J/ip + cc (dotted line) to J ftp momentum 
distribution for J/ijj production in T decay. And the sum of them is given by the solid line. 



On the experimental side, the CLEO collaboration fmd[24|] that the feed-down of Xcj to 
J/ip are < 8.2, 11, 10 percent for J = 0, 1, 2 respectively and the feed-down of ip(2S) is about 
24 percent in T — > J/ip + X. Therefore it indicates that the experimental result of direct 
J/ijj production would be 

£ Dircct (r -> J /ip + X) = 3.52 x 10" 4 (29) 

which is about 9 times larger than the presented theoretical results based on the color- 
singlet calculations. This means that unlike the conclusion before [26J the branching ratio of 
T — >• J/ip + X can not be explained by color-singlet model at the leading order. 

From the theoretical point of view, the color-octet mechanism can account for most 
J/ijj production, but its predictions for the J/ip momentum spectrum is not agree with 
the experimental data. The color-singlet predictions on the shape of the J/ip momentum 
spectrum is more closer to the experimental result, but the discrepancy of the branching ratio 
between them is large. For all the numerical results, we used the theoretically normalized 
decay width to estimate the branching ratio. Alternatively, by using (T|0 1 ( 3 S'i)|T) = 
2.9GeV 3 [2cJ to calculate the partial decay width and choosing the total decay width of 
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T 51.4 keV from the experimental measurement 34j], the branching ratio will be enhanced 
by a factor of about 3, which still can not explain the experimental results. Therefore, 
it means that the NLO QCD correction is important, just like in the known cases, the 
NLO QCD corrections for J /if) production in e + e~ annihilation show that the iCdactor 
are about 1.97 and 1.2 for e + e~ — > 7* — > J /if) + cc and e + e~ — > 7* — > J /if) + gg processes 
respectively; the NLO QCD correction in J/ if) related T exclusive decays are also found quite 



important 



351 ] . In addition, the contribution of 0(a®) processes bb( 3 Si, 1) — > cc( 3 Si, 1) +gg 



and bb( 3 S-\ , 1) — >• ce( 3 Si, 1) + gggg to the branching ratio has been estimated to be of 1CT 4 
order [25j. So that the next important step is to give an explicit and complete calculations 
of them, which will be very helpful to understand the conflict between the theory and 
experiment. Furthermore, to obtain the full QCD correction for the inclusive J /if) production 
in T decay would be a very interesting and challenge work for explaining the experimental 
data. But it will involve very complicated work at the QCD NLO and is beyond the scope 
of this work. 
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